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ABSTRACT 

We present results of combined N-body and three-dimensional reionization calculations to determine the 
relationship between reionization history and local environment in a volume 1 Gpc /i"' across and a resolution 
of about 1 Mpc. We resolve the formation of about 2 x 10* halos of mass greater than ^ IQ^^Mq at z = 0, 
allowing us to determine the relationship between halo mass and reionization epoch for galaxies and clusters. 
For our fiducial reionization model, in which reionization begins at z ~ 15 and ends by z 6, we find a strong 
bias for cluster-size halos to be in the regions which reionized first, at redshifts 10 < z < 15. Consequently, 
material in clusters was reionized within relatively small regions, on the order of a few Mpc, implying that all 
clusters in our calculation were reionized by their own progenitors. Milky Way mass halos were on average 
reionized later and by larger regions, with a distribution most similar to the global one, indicating that low 
mass halos are nearly uncorrected with reionization when only their mass is taken as a prior. On average, we 
find that most halos with mass less than lO'^M© were reionized internally, while almost all halos with masses 
greater than 1O'''M0 were reionized by their own progenitors. We briefly discuss the implications of this work 
in light of the "missing satellites" problem and how this new approach may be extended further 
Subject headings: cosmology: theory — galaxies: formation — galaxies: intergalactic medium 



I. INTRODUCTION 

The universe we observe at z = must bear the marks of 
reionization. Reionization began when the first stars pol- 
luted the intergalactic medium and created individu al H II 
region s (Alvarez et al. 2006; Abel et al. 2007; Yoshid aet al.l 
120071; IWise & A bel 2008). As the first galaxies grew in abun- 
dance, the H II regions became longer lived, eventually con- 
taining perhaps tens of thousands of dwarf galaxies, growing 
and merging until they overlapped, marking the end of reion- 
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Gnedin 2000a; Sokasian et al. 2001 


; Nakamoto et al. 2001; 


Ciardi et al 2003; Furlanetto et al. 


2004; lUev et al. 2006i; 


Zahn et al.1 2007t iTrac & CenI 12007) 


. Observations of high- 



redshift quasars imply tha t this proc ess was complete by red- 
shift z^ 6 (Becker 2001; lFanet al . 2002; W hite et al.ll200l 
IWillotl ? 2007), while large-angle polarization measurements 
of the cosmic microwave background constra in the duration 
of reionization (ISpergelll2003l: lKomatsull2008h . 

During this time, the temperature of the intergalactic 
medium increased from a few to tens of thousands of degrees, 
dramatically changing the evolution of gas as it responded to 
the highly dynamic underlying dark matter potential. 

Low mass halos in ionized regions are less able cool, col- 
lapse, and form stars than those in neutral regions, due to the 
increase in the cosmological Jeans mass when gas is ionized 
and photo-heated, someti mes called "J eans mass filtering" 
(e.g., [Shapiro et al. 1994; iTlioul & Wei nberg 1996; Gnedin 
l2000brTDiikstra et al.1 120041; IShapiro et a l. 2004). This sup- 
pression of structure is one of the fundamental ways that 
reionization can leave its imprint on subsequent structure for- 
mation, even up until the present day. 

Correlating reionization with the present-day environment 
may be the key to the so-called "missing satellite problem" — 
many more satellite halos are predi cted to form in CDM than 
are actually observed as galaxies dKlvpin et al.lll999h . The 
leading explanation — an alternative to more exotic possibil- 



ities like modifying dark matter or the amplitude of small- 
scale primordial density fluctuations — is that the UV back- 
ground maintains the intergalactic gas in a photo-heated state, 
preventing it from falling into the shallo w potential wells of 
the progenitors of th e satellite halos (e.g.. lBullock et al]|2000l: 
iBenson et al.ll2002l) . For example, one might expect that re- 
gions that were reionized earlier will have fewer luminous 
satellites than regions that were ionized later. However, biased 
regions, which are rich in early low-mass galaxy formation, 
would have reionized ^nsf. The latter effect implies that early 
reionization would lead to more satellite galaxies, while the 
former implies just the opposite. Detailed three-dimensional 
models are necessary in order to disentangle these compet- 
ing effects and quantify their dependence on the inevitable 
assumptions that must be made when modeling reionization 
on such large scales. 

In this Letter, we present our first calculations to ad- 
dress the correlation between reionization and local environ- 
ment. We take a novel approach, combini ng N-body simu- 
lations with a "semi-numeric al" algorithm dZahn et al.ll2007l: 
iMesinger & Furlanettoll2007h for calculating the reionization 
history for the simulation, allowing us to achieve a higher 
dynamic range in resolving the scales of reionization than 
has been possible until now. We then report on the statisti- 
cal correlations between halo properties and their reionization 
epoch and environment. We present our hybrid N-body /semi- 
numerical method in §2, our results in §3, and end with a 
discussion in §4. Throughout, we assume a flat universe with 
n„ = 0.25, CTg = 0.8, = 1, O,, = 0.04, and h = 0.7. 

2. MODEL 

Our hybrid approach consists of two steps. First, we run 
an N-body simulation of structure formation to determine the 
positions and masses of halos at z = 0. We then calculate the 
reionization history of the same volume in order to determine 
the reionization epoch of each halo. 
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Fig. 1 . — Visualization of the progress of reionization in our 1 Gpc/i"' calculation. Redshifts z = 14, z = 10, z = 8, and z = 6 are shown from right to left. Ionized 
regions are blue and translucent, ionization fronts are red and white, and neutral regions are dark and opaque. A random sampling of 5 per cent (about 40,000) of 
all the halos at z = are shown in yellow. Reionization is still quite inhomogeneous on these large scales, with large regions ionizing long before others. 



2.1. N-body Dark Matter Simulations 

For our cosmological N-body simulations, we used the code 
GADGET-2 (Springel 2005). We simulated a periodic box 1 
Gpc//; on a side with 1120^ N-body particles. We did not 
include any gas dynamics, a reasonable choice given that we 
are interested only in the global properties of the dark matter 
halos, not the internal properties of the baryonic component. 
We used a comoving softening length of 25 kpc, sufficient to 
resolve the formation of halos of mass ~ 10'^ M©. Atz = 
we use a friends-of-friends halo finder with a linking length 
of 0.2 mean inter-particle spacings to identify the halos. 

2.2. Semi-numerical Reionization 

Our mo del is based on the a nalytical formalism first intro- 
duced by Furlanetto et alj (12004,) and later ext ended to three 
dimensional realizations by IZahn et al.1 (l200l) . Its mam as- 
sumption is that a region is fully ionized if its collapse frac- 
tion, defined as the fraction of matter present in halos above 
some minimum mass M^in, is greater than some threshold. 
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This corresponds, for example, to the assumption that C./coii 
ionizing photons are released per atom per unit time. If re- 
combinations are neglected, then equation ([U results by en- 
suring that the time-integrated number of ionizing photons re- 
leased is greater than the number of the atoms. Another inter- 
pretation of the efficiency factor C, is that each halo produces a 
spherical ionized region around it, the size of which is directly 
proportional to its mass. Thus, all the recombination, and ra- 
diative transfer physics is absorbed into our choice of C,. For 
example, C, = ifescf*Nj/h)/i'i+nrec) where fesc is the escape 
fraction of ionizing photons from each halo, is the fraction 
of matter converted to stars within a halo, N^/t, is the number 
of ionizing photons produced in stars per hydrogen atom, n„'e 
is the average num ber of recombinations p er hydrogen atom 
during reionization (iFurlanetto et al.ll2004l) . 

To apply this criterion for "self-ionization" to an actual 
three-dimensional linear density field, we use the following 
relation for the collapsed fracti on within a sph erical region of 
mass m and density contrast S (iLacey & Colelil993i) : 
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v/2[a2(M,™„)-a2(m)] 
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where a^(m) is the mass variance over the scale m, Sdz) is the 
critical density for collapse, and M^i„ is the minimum mass of 



halos to be counted in the collapsed fraction, i.e. the minimum 
mass of a halo capable of producing a significant amount of 
photoionizing radiation. Note that the time dependence of the 
density field has been taken into account in the critical density 
for collapse, 5c(z) = Sc,o) /D{z), where D{z) is the linear growth 
factor, so that a{m) and S„ are constant in time. As shown 
by Furlanetto et al. (2004), this results in a time and scale- 
dependent "barrier" around each point. 
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6Am,z) = SAz)-V2[a\M^i,)-a\m)] ' err'(l-C-'). 

(3) 

The mean density within a sphere around a given point, 5m, 
must be greater than this barrier, 5_x{m,z), in order for the point 
at the center to be ionized by that region. A given point is 
considered to be ionized when the condition in equation [3] is 
met for any smoothing scale m, so that 



= MIN„ 
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where MIN,„ indicates the minimum value over all smoothing 
scales m. 

In practice, the outcome of this modeling is one value of 
Zreion at cach point on the grid, which characterizes the evolu- 
tion of reionization over all time. The smoothing of the den- 
sity field over all scales can be accomplished through a fast 
Fourier transform (FFT). We also store the radius at which 
each point on the grid first crossed the barrier, and associate 
it with it with the characteristic size of the region containing 
the sources that first ionized the grid point. 

Finally, to associate with each z = halo a reionization 
epoch and H II region size, we assign each of them a value 
that corresponds to the cell in which its center of mass lies at 
present. Given that typical H II regions are tens of Mpc, the 
vast majority of halos in our volume would not have had the 
required sustained peculiar velocities in excess of 10^ km/s 
for 10 Gyrs to have moved out of such a region. We thus ex- 
pect our results to be robust for most halos in the box, with the 
predictions being least accurate for the few halos that are just 
on the verge of falling into large galaxy clusters. Here, the 
reionization epochs could be overestimated, and H II region 
sizes underestimated. We set the parameters of the reioniza- 
tion model to have a minimum halo mass of M^ain = ^O^Mq 
and an efficiency parameter = 10. 

3. RESULTS 

Figure [T] shows snapshots of the reionization calculation at 
four different times for a 1024-' grid. Even on scales of 100 
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Fig. 2. — Reionization redshift distributions, for tlie regions tliat became 
halos in different mass intervals (see labels) by redshift zero. A uniform bin- 
ning was used in redshift, with the same spacing for each of the mass ranges. 
The low-mass distributions are very well sampled because of their large num- 
bers, while the high-mass halos, especially at the cluster scale, become noisy, 
owing to their small number (only 800 even in a 1 Gpc/i"' box). The thick 
red line shows the distribution for random points in the box. 

Mpc, reionization remains inhomogeneous, with the reioniza- 
tion redshift of regions as large as tens of Mpc varying be- 
tween Zr ^ ^5 and Zr ^ (i- Although regions that were reion- 
ized first are at the peaks of the underlying density field, there 
is not a one-to-one correspondence between the mass of the 
z = halos and their reionization epochs, since the halo and 
reionization barrier shapes and amplitudes differ. 

Figure |2] showns the distribution of halo reionization 
epochs, for several ranges of halo mass. There is consider- 
able spread in reionization redshifts in this model, ranging 
over 6 < z,- < 15. The most massive halos are biased toward 
higher values of z,-, peaking at z 10, 8, and 7 for masses 
M ~ 10'^ 10'^ and lO'^ M©. The distribution of the lowest 
mass halos, with masses ^ 10'^ Mq, does not have a well- 
defined peak, but rather increases toward the lowest redshifts, 
peaking at the percolation epoch at z ~ 6. This indicates that 
these lowest mass halos are relatively unbiased with respect 
to the structure of reionization. 

Figure|4]shows the 68 and 95 per cent contours of the reion- 
ization redshift distribution for halos binned by their mass. 
The median value increases from Zr ^ 8 for M/, = 10'^ M©, to 
z ~ 12 for Mf, = lO'^ Mq. The distributions have a long tail 
toward higher reionization values, which is more pronounced 
for higher-mass halos. Only 5 percent of lO'^-M© halos have 
Zr> 12, while only 5 percent of cluster scale halos have z,- < 8. 
Such a large spread in reionization epochs at all masses im- 
pUes that other halo properties, such as merger history and 
local matter density, may be important in setting the reioniza- 
tion epoch for a specific halo such as our own Milky Way. 

The mass-dependent distributions of bubble sizes are shown 
in Figure [3] Lower mass halos largely form in regions with 
larger H II bubbles, since their sizes increase with time. Inter- 
estingly, all of the roughly 800 cluster-mass halos in our sam- 
ple are associated with H II region sizes less than 30 Mpc. 
Only halos below about 10'^ Mq have H II regions sizes in 
excess of 100 Mpc, exceeding the mean free path for Lyman- 
limit systems and approaching (and potentially exceeding) the 



Fig. 3. — H II region sizes at reionization for the same mass bins as Fig- 
ure |2] The peak at low masses can be attributed to the grid resolution, given 
by a radius of ^ 0.5 Mpc. Increased resolution would reduce the effect, but 
on such small scales the role of recombinations and the stochasticity of the 
sources are likely to play an increasing role, effects we have not yet self- 
consistently included in these calculations. The distribution at scales larger 
than ~ 1 Mpc are robust predictions of our model. 

size of the box, 1 Gpc/i"' . 

4. DISCUSSION 

Using large-volume and high-resolution coupled simula- 
tions of reionization and halo formation, we have developed 
a new method for connecting the z = distribution of halos to 
the reionization epoch. We have found that, when only their 
mass is known, galaxy scale halos are nearly uncorrected 
with respect to reionization, with a distribution of H II region 
bubble sizes and reionization epochs that are roughly consis- 
tent with having a random spatial distribution. Higher mass 
halos, however, show a much stronger correlation, with none 
of the cluster scale objects having Zr < 8 or Run > 30 Mpc. 




10'= 10'= 10'* 



K^r [Mo] 

Fig. 4. — Median (solid line) and spread in the values of the reionization 
redshift, z^ion, as a function of halo mass. Shaded contours indicate the 68 
per cent (red) and 95 per cent (blue) spread in the distribution. The higher the 
halo mass today the earlier its progenitors were likely reionized. 
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EXTERNAL REIONIZATION INTERNAL FIEIONIZATION 

Fig. 5. — Schematic diagram indicating the difference between external 
and internal reionization. Ruii indicates the size of the reionizing bubble; 
Ring indicates the Lagrangian radius of the halo. 

An important distinction is between internal and external 
reionization. Figure |5] describes these two possibilities. In 
the external reionization case, the halo material was ionized 
by sources in a region with Rhu ^ ^Lag, where 7?Lag is the co- 
moving volume occupied by the mass of the halo at the cosmic 
mean density. In this case, most of the sources that ionized 
the material were not progenitors of the halo, and the ioniza- 
tion front swept over the halo's progenitors quickly, leaving 
the halo with a relatively uniform reionization epoch. For in- 
ternal reionization, Rhh « ^Lag, and the halo's reionization 
history is likely to be much more complex. In general, more 
massive halos were internally reionized, while less massive 
ones were externally reionized. 

Our defi nition is somewhat diff erent from previous defini- 
tions (e.g.. IWeinmann et al.l 120071) . but we believe our defi- 
nition is best suited for the method used here. In our defi- 
nition, halos are considered to be externally ionized if their 
Lagrangian radius, defined by Mhaio = 47r/y?Lag/3, is smaller 
than their H II region radius. For galaxy scale objects, with 
Lagrangian radii of order 2 Mpc, it is clear from Figure[3]that 
most of these objects were externally ionized. Because our 
model does not resolve scales below about a Mpc, however, 
it is difficult to determine how few galaxies were internally 
ionized. More detailed modeling of the small-scale structure 
on galactic scales, while still retaining the large volume pre- 
sented here, is therefore necessary. Our predictions for clus- 
ters are more robust. For these objects, with Lagrangian radii 
of order 20 Mpc, all but a tiny handful are internally ionized. 
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Our results may have important implications for galaxy for- 
mation, and in particular for the missing satellites problem. 
Because we find such a large spread in reionization epochs 
for Milky Way mass halos, more information, such as larger- 
scale environment and accretion history, will be necessary to 
determine the reionization epoch of our own halo, even after 
the global reionization history is well constrained. If the abun- 
dance of Galactic satellites is strongly dependent on the reion- 
ization epoch of the galaxy, our results indicate that Milky 
Way mass halos would have a large spread in the number of 
observable satellites. We explore this issue in a companion 
paper (Busha et al., in preparation). 

Our results may also have implications for the issue of 
galaxy "assembly bias", the idea that galaxy clustering may 
be dependent on properties other than the mass of the ir host 
halos (Wechsl eret al. 2006; Gao & White 2007; Crot on et all 
l2007h . If the reionization epoch of halos at a given mass is 
correlated with halo formation time, and if the reionization 
epoch affects any aspects of the galaxy population, then as- 
sembly bias could be more important for such galaxies than 
it is for their host halos. Further study will be required to 
investigate such effects. 

The approach we have presented here will serve as the foun- 
dation for such more detailed future studies. These studies 
will investigate the statistical correlations between present- 
day structure and reionization, and will also incorporate de- 
tailed galaxy formation modeling. Such improvements will 
allow for an investigation of the detailed coupling between 
star formation histories and the local reionization history. 
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